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Abstract 

The semiclassical theory for the large- iV field models is developed from an unusual point of 
view. Analogously to the procedure of the second quantization in quantum mechanics, the 
functional Schrodinger large-A equation is presented in a third-quantized form. The third- 
quantized creation and annihilation operators depend on the field V5(x). If the coefficient of the 
yj^-term is of order 1/N (this is a usual condition of applicability of the 1/A-expansion), one can 
rescale the third-quantized operators in such a way that their commutator will be small, while 
the Heisenberg equations will not contain large or small parameters. This means that classical 
equation of motion is an equation on the functional $[</)(■)]. This equation being a nonlinear 
analog of the functional Schrodinger equation for the one-field theory is investigated. The exact 
solutions are constructed and the renormalization problem is analysed. We also perform a 
quantization procedure about found classical solutions. The corresponding semiclassical theory 
is a theory of a variable number of fields. The developed third-quantized semiclassical approach 
is applied to the problem of finding the large-A spectrum. The results are compared with 
obtained by known methods. We show that not only known but also new energy levels can be 
found. 
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1 Introduction 



Investigation of quantum field theory as the number of fields tends to infinity is very important. For 
example, the number of colours of quarks Nc seems to be the only possible large parameter in quantum 
chromodynamics, so that one can perform asymptotic expansions in a series of 1/Nc |T], 01- 

Vector models are also suitable for the large-N analysis. Expansions in a series of 1/Nf {Nf - number 
of flavours of quarks) is widely used in QCD. Examples of applications of the 1/iV-expansion are: 

- large-N calculation of Green functions 0; 

- evaluation of the effective action 0]; 

- investigations of the spontaneous symmetry breaking 

- research of processes of particle creation and back reaction in the strong external field 0; 

- investigation of the evolution of the pion condensate 0; 

- evaluations of cross-sections of the processes like "1 particle" — > "n particles" 0. 

There are various approaches to the 1/A^-expansion. Consider, for example, the simplest vector model 

1 TD ^ A 

C = -d.v'^d.v^ - — y^V" - ^y^VVV (1) 

(a, 6 = 1, ...,iV, Lf^, are scalar fields, we sum over repeated indices). Some of these approaches are 

the following. 

1. One can extract the so-called bubble Feynman graphs and evaluate their contributions to the Green 
functions and scattering amplitudes P, 0|; one can then show that contibution of other graphs can be 
neglected. 

2. One can present physical quantities via functional integrals like 



J Dipe' 



fdxc^ (2) 



introduce the auxiliary field x Hj |T0|' HH ' insert the identity 

1 = J Dxe^ ^ 



into the functional integral (0) and evaluate the Gaussian integral over ip exactly. Then remaining integral 
over X is to be evaluated semiclassically (by the saddle-point technique). It is the condition N ^ oo that 
allows us to apply the semiclassical approach. 

3. One can use the collective field approach |12[. The Schrodinger wave functional \E'[(y9i(-), (y9jv(-)] 



can be assumed to depend only on the collective field 

N 

^(x,y) = ^^'^(x)^'^(y). (3) 

a=l 

The ansatz 

vl>[y,i(.),...,¥.^(-)] = ^[^(-,-)] (4) 

satisfies the Schrodinger equation; the obtained relation for \E'[V^(-, ■)] can br treated semiclassically. 
4. One can develop the following semiclassical approach [|T3]. Consider the operators 

N 



7A(x,y) = 5:r(x)(^'^(y) 

a=l 

N 

ax,y) = Enx)r(y) 



a=l 



1 



N 

r)(x,y) = ^7r'^(x)7r"(y). 

a=l 

The system of Heisenberg equations for these operators is expressed via these operators only, so that one 
can use the argumentation hke the Ehrenfest theorem in quantum mechanics. One can consider the average 
values of the right- and left-hand sides of the Heisenberg equations and obtain the classical equations. 

5. One can consider the Heisenberg equations for the fields (f"'. Making use of the averaging procedure 
[ill, one can obtain the large-N equations. 

In this paper we present yet another approach to the 1/A^-expansion. We show that semiclassical 
investigation of the large-N theories may lead to classical equations of another type rather than to equations 
on "collective field" ip{'x,y) or yj'*. 

We consider states of a more general form than "^[ipl-, ■)]. The state functional 

vl/[^i(-),...,^^(-)] (5) 

is considered to be symmetric with respect to transformations of fields v'i(x), ip]^{x.). We will also 
discuss the general, non-symmetric, case. 

In quantum statistical mechanics the most suitable way to investigate the iV-particle Schrodinger equa- 
tion is the method of second quantization. The Hamiltonian operator is presented via the operators of 
creating and annihilating particles, a^{q), with given coordinate q. Making use of the canonical commu- 
tation relations, one can investigate different properties of the A^-particle system. 

One can hope that this conception can be applied also to the iV-field system. The A^-field Hamiltonian 
operator is to be presented via the operators [</)(■)] which "create" or "annihilate" the field configuration 
This procedure can be called as "third quantization" because quantum field is already second- 
quantized. 

The notion of third quantization usually arise in quantum gravity and cosmology [|T^ when the processes 
(like wormhole transition) with variable number of universes are studied. Creation (annihilation) operators 
create (annihilate) the universe. 

The conception of third quantization is also useful in string theory: splitting and joining of strings are 
interpretted in terms of creation and annihilation operators ||T6[| . 

We see that the idea of third quantization can be applied also to ordinary field theories. We will show 
that this idea is very useful and allows us to construct new asymptotics in large-N field theories. 

For example, one can investigate the processes with a large number n of particles, n ~ A^, by using the 
developed technique. It is known p that the usual 1/A^-expansion is applicable only for processes with a 
small number of particles. 

We will consider the formalism of third quantization in more details in section 2. We will show that the 
coefficient of ip^ in eq.(|l|) is of order 0{1/N) is a condition of applicability of semiclassical approximation 
to the third-quantized theory. The corresponding classical equations are unusual. Namely, the classical 
equations corresponding to quantum mechanics are ordinary differential equations for the classical trajec- 
tory {p{t),q{t)), p{t) are momenta, q{t) are coordinates. When one applies semiclassical approximation 
to the (second-quantized) field theory, one finds classical equations for the classical field ip{'x,t). In the 
third-quantized theory the classical equation is equation for the functional ip{-)] depending on the field 
(/9(x) and time t. Although this equation resembles the ordinary quantum functional Schrodinger equation 
in the theory of one field, an essential feature of this equation is non-linearity. Another interesting feature 
of the equation on </)(■)] is exact solvability for some interesting cases. 

The conception of classical "master field" is often used in investigation of large-N systems. However, 
"master field" is usually considered as a function on a finite-dimensional space. For example, field '?/'(x, y) 
(H) or fields 93" (x) are ususally considered as classical variables [|l^, |l^]. We see that in our approach 



the role of "master field" is played by the functional </?(■)] which can be treated as time-dependent 
vector of the one-field Fock space. Classical mechanics for '?/'(x, y) is in agreement with our dynamics for 
$. Namely, the Gaussian ansatz for our equation for $ leads to the equation on ip which was obtained 



2 



in refs.|jT2|. However, our classical theory is much more rich than the mechanics for ip because one can 



consider other solutions for the equation for $. The reason is that we are considering solutions to the iV- field 
equation which are not 0(iV)-symmetric but symmetric with respect to transpositions. This class is much 
more wide than the class of solutions This means that we are able not only to reproduce the known 
asymptotics but also to construct new asymptotic solutions to the large-N quantum-field Schrodinger 
equation. 

An interesting feature of the classical equations on is that it can be interpretted as a Hamiltonian 



system corresponding to the flat phase space. On the other hand, the phase space of refs.|13] is curved 



There are a lot of well-developed semialassical methods |T^, |2^, ^ for flat phase space, so that all 

of them are applicable for the third-quantized system. 

We will construct the asymptotics for the entire N- field wave functional (^. This approach gives 
more information than the averaging procedure of the Heisenberg field equations. The difference can be 
explained as follows. Consider the semiclassical approximation to the ordinary quantum mechanics. One 
can derive classical equations with the help of the Ehrenfest theorem. This result allows us to find the 
trajectory of propagation of the wave packet. To find the evolution of the shape of the wave packet. 



it is necessary to apply more complicated semiclassical methods such as the complex- WKB method |^ 
which allows us to construct the approximation for the wave function rather than for average values of the 
semiaclassical observables. 

This paper is organized as follows. Section 2 deals with the simple quantum mechanical model, the 
0(A^)-symmetric anharmonic oscillator. One can test different versions of the 1/A^-expansion on such a 
toy model. We consider the analogs of the collective field approach and of our approach. These methods 
are compared, and it is shown that our approach allows us to construct much more rich set of approximate 
solutions to the Schrodinger equation. In section 3 we introduce the notion of third quantization (as 
related to the large- system). We write down the third-quantized quantum Hamiltonian, as well as 
the Heisenberg equations of motion. The semiclassical approximation is applied and the problems of 
regularization and renormalization are considered. In section 4 we apply the developed method to the 
problem of investigation of the large-A^ field theory in the finite volume. Section 5 contains concluding 
remarks. Appendices A and B deal with the brief review of the complex- WKB technique 3S 



2 0{N) - symmetric anharmonic oscillator: different approaches 
to 1/A^-expansion 

In this section we consider the simplest example of the large-N system - the 0(A^)-symmetric anharmonic 
oscilator with quartic interaction. The wave function of this system depends on A^ coordinates, xi,...,xn, 
while the Hamiltonian is 

Different approaches to 1/A^-expansion can be illustrated by applying to such a simple system. Let us 
compare the third-quantized approach developed in this paper with other approaches. 

2.1 The collective field approach 

The collective field approach for the Hamiltonian which is analogous to the ansatz (D was developed 



m 



r2| . The idea is to consider the 0(A^)-symmetric wave function depending only on r = -i/xf + ... + 



Af5 
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However, the full probability will not have the usual form / dr\ilj{r)\'^, since the area of the surface of the 
sphere depends on its radius. This means that 



dxi...dxN\i'N{xi,...,XN)\'^ = Cn c/r|V^(r)|V^ \ 
where Cat is an r- independent constant. Let us consider the quantity 

N—l 

Lp(r) = il){r)r 2 , (7) 



which plays the role of the probability amplitude that r = yxf + ... + x^. Time evolution of the function 
(1^ y9*(r) is specified by the equation derivable from eq.(|D: 

.dip 1(9^09 A^^ - 1 ar^ q . 

I— = (pH 09 + — rV (8) 

dt 2 ^ 8r2 ^ 2 ^ 4Ar ^ ^ ' 

The rescaling 

r = ^^/N 

transforms eq.(^ to the semiclassical form 

+ ^.{^-T7^.)v + ^'P + -rV^ (9) 



N dt 2A^2 5^2 8^2^ p^2j^ 2^ A 

where the analog of the Planck constant is 1/N. The condition — 00 is the condition of applicability 
of semiclassical methods to eq.(|^). 

Eq.® specifies the motion of the particle in the potential 

One can apply different semiclassical methods to eq.(^ which are described in Appendix A. Spectrum in 
the vicinity of the ground state can be found by using the oscillator approximation of refs. : 

^ = UiO + l^^JuW){n + 1/2) + 0{1/N'), (11) 
where ^ is the minimum of the potential (p!0|). 

2.2 The "second-quantized" approach 

Let us develop for the 0{N) - symmetric system the approach based on the Fock space. Instead of the 
Schrodinger equation for ip^^xi, ...,xn), we consider the evolution equation for the set of functions 

(V'o, V'i(a;i), ...,^pN{xl, ...,xn), ■■■) (12) 

{ipN being symmetric). Introduce in the Fock space of sets (0) creation and annihilation operators as 
usual, 

(a"(a;)^)„_i(xi, ...,x„_i) = ^/n'^|Jn{xl, ...,a;„). 



1 " 

(a+(x)V^)„(Xi, X„) = ^ V 6{X - Xi)V'n-l(Xl, Xi-i, Xi+i, Xn). 



Since the operator 



E 

1=1 
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is presented via the introduced operators as 



j dxa'^{x)x'^a (x), 



while for the operator 9^/i9x^) one has 



dxa~^(x) 



dx"^ 



a (x) 



one can present the Schrodinger equation in the second-quantized form 



where 



H 



dxa^{x) 



1 9^ ax^ 
'2dx^ ^ ~Y 



a (x) + ^ [ I dxa^ {x)x^ a {x) 



(13) 



(14) 



for e = 1/N. Since the operators a^(x) obey canonical commutation relations 

[a^{x),a^{y)] = 0, [a"(a;), a+(y)] = 6{x - y), 

one can use the following functional Schrodinger representation. The vacuum state of the Fock space 
(1,0,0,...) corresponds to the functional 



^^J{^) = exp(^-^|dxeix: 



while creation and annihilation operators are written as 

ax)Te6/6ax) 



a (x) 



2e 



In this representation eq. (|13]) takes the semiclassical form 

6 



where 



Ot 



5^ 

^ — in ^ + m 



(15) 



(16) 



E{ip\^\ = J dxip*{x) (^-l^ + ^) + !(/ dxx^\ip{x)\'- 

An analog of the Planck constant is e. One can try to apply semiclassical methods to eq. (p!5D . One finds 
that the corresponding classical equations 



d_^_5H d_^_ _6n 
dt 6tt ' dt 6^ 



(17) 



can be presented as 



d 



SH 



dt Sip* 
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We see that classical trajectory in our approach is specified not by one coordinate and one momentum but 
by infinite set of coordinates and momenta, (,^(-), 7r(-)), or, equivalently, by the complex function (p{t,x) 
obeying the following dynamical equation 

= (-^^ + "-^Mt,x) + / dyy^\^{t,y)\\ (19) 

The reason for extending the phase space is the following. First, we consider not the "N- argument" 
Schrodinger equation but infinite sets of such equations for = 0, oo. The average value of "number of 
arguments of wave function" analogously to the average number of particles in quantum statistics 



< J dxa^{x)a {x) >= ^ ^ J '^^ ^^(^) ~ g ^^^ ^ (^^i^) + ^-rr^ 1 > 



in the "semiclassical state" is approximately equal to 



^ J dx{^{x) - in{x)){^{x) + iTT{x)) = ^ J dx\i^{x)\'^. 
Because of the relation e = 1/N one should impose on the function (/? the normalization condition 

dx\(p{x)\'^ = 1. (20) 



However, the extension of the state space to the Fock space (|T2]) is not the only reason for extension 
of the phase space. Another reason is that we consider not only 0{N) - symmetric solutions to the 
Schrodinger equation for ijjpf but solutions of a more general form. 



2.3 Correspondence between Hamiltonian systems 

Let us show that the collective-field equations can be obtained as a partial case of our classical equation 
(|IDp. To investigate this problem, consider the substitution 



y,(t,a;) =ce^"^' (21) 

to eq.(|19D for t-dependent complex numbers c and a. Making use of the condition (^), we find that the 
function (^ obeys eq. (|l^) if 

'dt'=-^^ (2'^ 

^ + «^ + «+27b-° (23) 
However, eq. (|23| ) can be transformed to the form of the Hamiltonian system corresponding to the potential 
(|TUp. Let us extract real and imaginary parts of a, 

Rea 1 



\pTrma ^/21'ma 
Eq. (|23|) takes the form 

dx dp dU 

dt ^' dt dx 

We see that our approach contains the collective field approach as a partial case. However, one can consider 
not only the Gaussian ansatz (21) for eq.(19) but also other substitutions. We see that one can construct 
new asymptotic solutions for the large-N systems. 
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2.4 Asymptotic spectrum of the second-quantized Hamiltonian 

Let us show that the oscillator approximation applied to the Hamiltonoan system with the Hamiltonian 
(|16D leads us to the semiclassical spectrum of energy which is much more rich than eq.(|lT]). 

There are many semiclassical methods. One of them - quantization of periodic trajectories - is the 



following [21|. One should find the classical periodic solution. Then one considers small variations around 



it and finds the "stability angles" |2T| related with the distances between semiclassical energy levels. This 
approach is certainly applicable to the Hamiltonian system (|I6|) since there are many periodic solutions 
to eq. (|T9D of the form 

(p{t,x) = ip{x)e-'^'. (24) 
However, we will apply the much more simple approach - quantization of static solutions. To apply this 



approach, it is necessary to consider the Hamiltonian 

Hfi = H — n I dxa'^(x)a~ (x) 



(25) 

instead of the operator (0). Since the spectra of the A^-particle parts of these Hamiltonians are related 
by shifting by QN, one can investigate the Hamiltonian (|25|) instead of (|T4|). 

The classical Hamiltonian function corresponding to the operator (^) has the form (p!6D, where 



Hn[ip*,ip] = J dx(p*{x) 



1 d'^ ax"^ 
'2dx^ ^ ~Y 



n\ Lp{x) + I dxx'^\ip{x)\' 



(26) 



Let us find the .static solutions of the classical equation of motion (|T8D which takes the form of the 
stationary harminic oscillator equation 



1 



oj'^x'^ 



Vl (p{x) 



\ 2dx'^ 2 

where the frequency u is expressed via the average value of x^: 

a + g J dxx'^\(f{x)\" 



0, 



(27) 



f28) 



Note that the same equation ([27|) could be obtained by substitution of periodic solution ( ^4]) to the equation 
The solution to eq. (pT)) can be chosen as the i^-th excited state of the harmonic oscillator: 



K 



9? = \E' 



K 



xfK'\ 



0; 



where 

is the ground oscillator state, while 



A' 



conste 2 

1 d 



X 



2uj dx 



is the creation operator. The parameter Q is 



The normalizing factor is determined by eq 
creation and annohilation operators. 



Q = uj{K+l/2). 

I). Since the operator of coordinate can be presented via 
1 



X 



{A+ + A- 



(29) 



one can find the average value 
so that the equation on oj (^) can be presented as 



cu^ = a+-iK + 1/2). 

Let us quantize the constructed static solution. It is necessary to consider the small perturbations 
around it: 

^(2^) ^ ^(3;) + ^^(2^? T^ix) 7r(x) + 57r(x, t). 

Then one should consider such variations that 

t) = 6^{x)e'^\ 67r{x, t) = 6n{x)e'^'. 

This means that variations of classical variables ^ and vr should be complex, although the classical variables 
are real. It is convenient to make a transformation: 

The fact of complexity 6^ and 6tt means that the variations of and (p* are not conjugated each other; 
they are independent. Thus, to investigate the small variations around static solution, one should consider 
the Hamiltonian system 

.dip^6Hn .dip* _ 6Hn 

dt 5(p* ' dt 6(p '' 

consider the independent perturbations of ip and ip*: 



and find a spectrum of /5: 



+ Fe'f"', ^* -^'p* + Ge 



6ip*6ip 6ip*6ip* ' 



6ip6ip 6ip6ip* 



(30) 



which should be real if the static solution is stable. The more detailed derivation of the semiclassical 
spectrum is presented in Appendix B. 

The variation system (^) has the following form for the Hamiltonian (^): 

-W^)= + Fix) + ^-x'if'ix)JdyyMy)iF + G)iy), 

(31) 

= (-1-^ + ^ - ^) ^(^) + f ^V'(x) / dyyMy)iF + G){y), 



The solutions to eq.(|3lD are presented as follows. First of all, let us try to substitute instead of F and G 
the n-th eigenfunction of the oscillator: 

F+ = 0,a+ = ^„, (32) 
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or 

F- = ^„,G';=0. (33) 
Since the matrix element {'^k, x'^'^n) = if n — K ^ 0, ±2, functions (|32|) , (|33D obeys eg .(^11) if 



Let us find other solutions to eg. (|3lD . First of all, there is a "zero- mode" solution corresponding to the 
invariance of the Hamiltonian system with respect to transformations if — > 9?e*", if* — >■ 99*6"*": 



Further, making use of eg. (p9D , one has 



Thus, we find the following two solutions to system (PT|): 



F = JK{K - l)^i,+2, G = -J{K + 1){K + 2)^K^2, (3 = -2a;; 



G = ^K{K - l)^K+2, F = -^{K+1){K + 2)^K-2, f3 = 2uj 

because for these cases {x'^'^k, F + G) = 0. 
The non-trivial solution to eg.(|HT|) is 



^ - ^^^^ " (3-2u 



^'{K+l){K + 2) ^K{K-l) 
where 

= Au^ + ^{2K + 1). 

There is also a solution to the non-stationary variation system which is proportional to t: 

5^* -b^ = 2iuj{2K + l)t^K, 

3 



S^* + 5ip = -^iK + l){K + 2)<ilK+2 + ^K{K -l)mK-2 + (2 + J^RTV)^^'^ ' 

However, such linear instability of classical static solution is usual for guantum systems with zero-modes. 
There are several approaches to resolve the difficulty. It is shown in Appendix B that such an instability 
is not an obstacle for guantization. 

The asymptotic energy spectrum for the guantum Hamiltonian Hq^ is expressed via freguences j3: 

= NSn + 60 + + 0{1/N), (34) 

m 

where is a guantum energy, is a corresponding "classical" energy, 

Eq is a guantum correction, are energy excitations. Non-negative numbers Um are parameters specifying 
the energy levels. 
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It is important that not all values of /3m found above enters to eq- fplf) . When we quantize ordinary 
harmonic oscillator (or free quantum field), there are positive-frequency and negative- frequency solutions 
to the variation equations (for free fields they coincide with classical equations). However, one usually 
takes into account only positive frequences to construct a spectrum. 

It happens (see appendix B) that in our case we should take into account not positive values of P but 
such values of f3 that obey the following condition: 

(G,G)>(F,F). 

This means that the excitation energies entering to eq.(p^ are: 

pm = ujim - K),m^ K,K ±2, 

pK-2 = -2cu, 



PK+2 = +xku;^ + -{2K+l). 
V uj 

An asymptotic energy spectrum of the operator H which corresponds to the pertiodic solution ( p^ ) has 
then the form E = Eq + NQ, or: 

E^ = NS + eo+ E Pmnm + 0{1/N), (35) 

m>0,my^K 



where classical energy is 



S = H{^\ = ooiK + 1/2) -^,{K + 1/2)1 



2.5 Correspondence between asymptotic spectra 

Let us investigate the correspondence between found asymptotic formulas for the spectrum, eqs.(|l^) and 
(pSj). Let us set in eq.(p5D K = 0. In this case the frequency u obeys the equation 

uj^ = a + ^. (36) 
Let us find the minimum of the potenrial (plQl). The corresponding equation for ^ is 

-^ + a^ + gt = 0. (37) 

However, eqs.(0) and (|37|) coincides if one sets: 

f = ^. 

^ 2u} 



It is easy to see that £ = U{^), while = yU"{^). Therefore, the asymptotic spectrum in the case K = 
is 



= NUiO +^0 + \/U"iOn2 + ujni + Sljus + ... 
The result can be interpretted as follows. There are many perturbations around classical solution. However, 
the excitations with energy do not break the quantum 0(A^)-symmetry, while other excitations break 
this symmetry. This is a reason that they have not been discovered by the collective-field approach which 
allows us to find 0(iV)-symmetric solutions only. The method of second quantization allows us also to 
find 0(A^)-asymmetric solutions. 

We can also notice that one can perform a quantization around solutions with K ^ 0. These solutions 
do not have form (^Tj), so that the corresponding quantum states are 0(A^)-nonsymmetric. One can say 
that quantization of the K = 0-solution is investigation of symmetric and nonsymmetric perturbations 
around 0(iV)-symmetric quantum state while quantization around other solutions correspond to research 
of the fluctuations around non-symmetric state. Thus, the second-quantized approach allows us to find 
new energy levels of the anharmonic oscillator. 
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3 Formalism of third quantization and semiclassical approxi- 
mation 



This section deals with the third-quantized formulation of the quantum- field theory ([I|). We present the 
Hamiltonian operator of this theory via the creation and annihilation operators A^{(f{-)). We are also 
going to consider the problem of applicability of the semiclassical conception to the third-quantized theory. 
The "classical" equations for <l>[(y9(-)] are to be derived. 



3.1 Creation and annihilation operators 

In the functional Schrodinger representation states of the N-field system are specified by (time-dependent) 
functionals (W), while their evolution is described by the functional Schrodinger equation 



^—7— = Hn'^n, 
at 



(3J 



where 



N 



(ix 



1 52 1 \ 



+- 



dxv^a (x) (x) (x) V5b (x) , 



where e is a coupling constant in the A^-field theory, e = 1/N. 

The first step to construct the third-quantized formulation is to extend the Hilbert state space. Instead 
of the space of functionals \i/[(y9i(-), (^7v(')]) l^t us consider the "extended state space" which is a Fock 
space of sets 

/ *o \ 



*iv[<^i(-), 



•)] 



\ 



(39) 



where is symmetric with respect to transpositions of ipi and (pj. One could conclude that the physical 
meaning of \&0; ^ivj^Afv is the following: ^0 is the probability amplitude that there are no fields, \E'i 
is the amplitude that there is only one field etc. However, theories with uncertain set of matter fields are 
not of widely use, so that one should consider only such extended state vectors (|39|) that ^ and 
\E'fc = 0, k ^ N . However, considering the extended Schrodinger equation 



' dt 



Hk'^k,k = 0,1,... 



(40) 



on the Fock vector ( [59| ) instead of eq. ( pq) will be shown to be very useful, since system of equations 
is more suitable for the semiclassical analysis. One can construct asymptotic solution to the system 
and then extract the A^-th component of the extended state (P^l). 

Eq.(POD can be presented via the creation and annihilation operators. These operators A^{ip{-)) trans- 
form the extended state (|39|) to the Fock vector 



(A^ivi-WhiM-)) 



\ 



{A^ivi-mMi^i 



.<^7v(-)) 



V 



(41) 
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where 



1 ^ 



V a=l 



(42) 

The functional 5-function entering to eq. (^) can be formally defined via the functional integral over all 
spatial field configurations 

jD<p6{<P{-)-y.{-))F[<p{.)]=F[v{.)] 
for arbitrary functional F. The formal integration measure enters also to the inner product, 



The introduced operators (^) obey usual canonical commutation relations 

[A±[v.(O],A±[0(O]] = O,[A-[^(-)],A+[0(-)]] = %(-)-0(-)) (43) 



Let us present different operators in the extended state space ( p9D through the operators (|4^) . It follows 
from eq. ( ^2]) that the operator of "number of fields" which multiplies the k-th component of ( pOj) by k can 
be presented as 

' D^A+[^{-)]A-[^{-)]. (44) 



Analogously, the operator 

k 

a=l 

can be written as 

J DifA+[ipi-)M^My)A-[cpi.)]. (45) 

One has also the following relation 

E , . , = I D^A+[^{.)]--^-—A-[^{.)]. (46) 

S'^aip^maiy) J (5(/?(x)5v?(y) 

Making use of the obtained relations, we present the quantum Hamiltonian via the creation and annihila- 
tion operators 

jD^D4> I dx^'^(x),#,2(x).4+b(.)lA-b(.)l-4+|,#.(.)l-4-W01. (47) 
The system of equations (^) can be presented in a simpler form 

= H^, (48) 

at 

where \E' is an extended vector (pOj), H is the operator (|47|). 
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3.2 Perturbation theory or semiclassical approximation? 

Let us consider the case of small coupling constant e. At first sight, at small e the theory ( ^TD can be 
approximated by the theory of free fields. One can use then the perturbation theory which can be shown 
to be equivalent to the ordinary, second-quantized, perturbation theory. 

However, we know from the quantum field theory that there are non-perturbatve effects even at small 
values of coupling constant. Namely, the conditions of applicability of such non-perturbative methods 
as soliton quantization pll], quantization in the background of the strong external field EH] or instanton 



method |2^, ^ is that the coupling constant should be small. 

To investigate the problem of applicability of perturbative or non-perturbative approach, consider the 
Heisenberg equations for the operators 



which have the form 



dt' 

and can be simplified by using the canonical commutation relations (^3l). The equation for the annihilation 
operator is 

d /■ / 1 (5^ 1 \ 

+y Jd^, I dx^Hxmx)Aim]Aim]A;if(-)] m 

the conjugated equation is the equation for the creation operator. 

The usual procedure to derive the classical equations of motion is averaging of eg. (^91) . If the second 
term of the right-hand side of eq.(|i9|) can be neglected, the perturbation theory is applicable. Otherwise, 
one should develop a non-perturbative approach. The difficulty may arise from the factor 

|d0A+[0(-)]0^(x)A-[0(-)] 

which is of order k when it is applied to the fc-field component. Thus, we see that if the number of fields 
satisfies the condition 

« 1/e 

the perturbation theory is applicable, while if 

A^ ~ 1/e, 

another approach is necessary. 

3.3 Classical equations 

3.3.1 Averaging the Heisenberg equations 

One of the approaches to the semiclassical theory is the following. Let us rescale the creation and annihi- 
lation operators 

= V^A^M-)]- 

The Heisenberg equation (l49|) will take the form 



i-t-M-)) = y [-^j^f^^^ + 5V^(x)v^(x) + -^(x)^(x)j 
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+ ^Id<pI dx^2(x)02(x)$+(</.(-))$,-(</>(-))$r(v^(-)) (50) 

which does not contain the small parameter e. 

The coupling constant e appears, however, in the canonical commutation relations 

[^tH-)],^t[m] = em-)-v{-))- (51) 

The fact that the small parameter arises in the commutatora and does not arise in the dynamical equa- 
tions resembles semiclassical quantum mechanics: the Heisenberg equations are regular in semiclassical 
limit (they transforms to Hamiltonian system), while the commutator between coordinate and momenta 
operator tends to zero. 

Analogously, one can try to substitute the operators entering to eq. ([50|) by the c- number quantities. 



and obtain the classical equation 

d /■ / 1 5^ 1 \ 



+^ Id4, I dx,p2(x),#,2(x)4.,-(,#,(.))*,W0)*,(43(0) (52) 



Such substitution can be explained via the averaging procedure of eq. ( |50| ) analogously to the derivation of 
the Ehrenfest theorem in quantum mechanics. Since the commutator between operators is small, the 
uncertainty of them can be made of order 0{y/e). This means that for some states average value 

<Ai^t,^;)> (53) 

can be approximated up to 0{y/e) by the quantity A($f , ^t), 

<Ai^t,^;)>^Ai^;,^,). (54) 

This means that classical limit of eq . (^OD is eq . (^21) . 

The states obeying eq. (^4]) can be constructed as follows. Let Y be e-independent extended state (p^). 
Consider the unitary operator 



f/$ = exp 

obeying the relations 



/ D<pmv{.))A-^{v{.)) - <i.*(^(-))A-(^(-))) 



(55) 



U^^A+U^ = A+ + ^^^^ 
which are corollaries of the canonical commutation relations (^SD. Consider the semiclassical state 

U^Y. (57) 



It follows from eqs. (|56D that average value (|53|) is equal to the average value of the operator A{^^ + 
^At,^t + V^A I ) over e-independent state Y . Making a limit £ — > 0, we obtain eq.(B3). 

Thus, we see that semiclassical analysis of the A^-field equation based on the Heisenberg approach may 
lead us to classical equation (0). The classical variable is a complex fiunctional $[</)(■)] being the state of 
quantum one-field configuration. This state obeys unusual nonlinear equation. 
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The "semiclassical" state (|57|) is a Fock vector ( pOf ) with all non-zero components. One can try to 



extract from eq.(|5^) one of the component, for example, the A^-th component. However, one should be 



careful because our calculations are not exact but approximate. If the A^-th component of eg . (|57|) is 
exponentially small as e — 0,it should be neglected, so that the approximate wave functional will vanish, 
and the constructed asymptotics will be trivial. For the simplest case {Y is a vacuum) the probability 
that there are N fields is given by the Poisson distribution with the maximum at 

N = yDm<pm'- (58) 

This relation can be derived also from eqs.(Q) and (0). Eq.(|58D means that semiclassical approximation 
allows us to investigate the large-N field theory if ~ Note that the case iV ~ 1 corresponds to 
$ ~ ^/e, classical equation (^) transforms to the equation for free fields, so that the perturbation theory 
is applicable. 

3.3.2 BBGKY-like approach 

Eg . (|5^) can be also derived by the BBGKY (Bogoliubov - Born - Green - Kirkwood - Yvon) approach 
(p8|], see also |^) from the A^-field equation (^). Let us briefly discuss this technique. 
Consider the fc-field correlation functional 



Rki^ii-), ^k{-); 0i(-)> M-)) = J Dipk+i...D(pN 



X^'^[y^l(-), ipNi^NlM-), M-), Vk+li-), ^Ni-)] (59) 

corresponding to the A^-field wave functional "^n. The correlator (|59[) can be presented via the creation 
and annihilation operators 

{N-k)\ 



< A^[M-)]...A^[M-)]A-[M-)]-A-[M-)] > 



It follows from eq.(49) that functionals (|59D obey the BBGKY-like hierarchy of equations. 



Suppose that the correlators ( p^ factorize as k = const, N ^ oo 

Ri ~ $:[^i(-)]...$n^fc(-)]'^'t[0i(-)]-<^'*[0i(-)]- (60) 

Substituting eq.(0) to the BBGKY-hierarchy, one can find that for some phase factor e*'''* the functional 
^te'^' obeys eq.(||). 

3.3.3 Variational principles 

The classical equation of motion ([5^ ) can be presented as a Hamiltonian system. Namely, this equation is 
derivable from the variational principle 



t 

dt 







D.pfM-)]ijtM-)] 



2 (5(/9(x)5(/)(x) 

^ I D^D<P I d^^'i^)^'i^)\<!>M-Wmm\' ] - extr. (61) 



This variational principle is Hamiltonian. Namely, one can extract real and imaginary part 

$ = — ^— 
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and notice that the principle ( |6ll) is 

J dt{PQ - H) ^ extr. 

Thus, classical equations of motion correspond to infinite-dimensional fiat phase space rather than to 
curved space considered in ref . [l^ . 

We will see that the (time-dependent) Gaussian ansatz 

<l>[<^(-)] = cexp Q j (^xdy<^(x)A(x, yMy))) (62) 

satisfies eq.(0), so that one can consider the subspace of the phase space which is associated with Gaussian 
states (p^). The corresponding variational principle for classical mechanics can be obtained by substitution 
eq. (|6^) to eq. (|6T|) . The obtained dynamics coincides with refs. ||T3|, |T^. 

However, one can consider other solutions to eq.(^), for example, the product of the polynomial in 
by the Gaussian exponent (|62D. These solutions will not belong to the reduced curved phase space but 
will belong to the fiat space {$[(y9(-)]}. 

The variational principle (|6TD can be interpretted as follows. One can consider the time-dependent 
variational principle ^ 

J^dt{^%,{zj^-HN)^%) (63) 
for the A^-field equation and substitute the A^-particle time-dependent test functional 

vl>*^[¥,i(-),...,¥^^(-)] = $tbi(-)]-'^'*b^(-)] (64) 

to eq. (|63D . One obtains then the classical variational principle (|6T1) . 

Note that the Gaussian ansatz to the quantum variational principle instead of 
mechanics in the reducrd curved space. 

Although the wave functional (|6^ will be shown to be not the asymptotic solution to eq . (|38D , eq.(^) 
is correct and can be derived also by the substituting the asymptotics for the wave functional to eq. (p8|) . 
This will be done in section 4. 



34) lead to the classical 



3.3.4 Operator formulation of classical mechanics 

Let us discuss relation between eq.(^2|) and Heisenberg approach ||l4l. Consider the "Heisenberg" field 
(^*(x) which can be defined as follows. Let be a non-linear operator transfering initial condition for 
eq.(^) to the solution to eg. (|52D at time t: 



The Heisenberg operator is 
such that 

It follows from eqs.(p2D and 



^*(x 



($,^*(x)<l>) = ($*,<^(x)$*). 
that the Heisenberg field <^*(x) = ip{x) obeys the equation 



d^d^(p{x) + w?(p{x) + A < (^^(x) > (p{x) 







(65) 



(66) 



which was obtained in 14 



An interesting feature of this second-quatized theory which is classical for the A^-field model is that 
the field ^{x) is a non-linear operator in the one-field state space. We will discuss in the next sections 
the applications of the semiclassical approach to eq.(p8D and clarify the role of eq.(^) in constructing 
asymptotic solutions to the A^-field equation. 
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3.4 Regularization and renormalization 



In the previous subsection we have derived the classical equation ( |5^ on the functional However, 
this derivation was formal. We have not taken into account the problem of divergences and renormalization 
in quantum field theory. Thus, eq.(E^) is not well-defined, and it is necessary to investigate the problem 
of correct definition of the classical equations in details. 



3.4.1 Regularization 

There are many ways to regularize the quantum field theory. Relativistic-invariant regularizations are 
usually applied to evaluation of Feynman graphs, while lattice regularization is studied in non-perturbative 
approaches. Another way to regularize the theory is to substitute the field V5(x) by the cutoffed field v5a(x): 

(/?a(x) = J c/yAA(x - y)^(y), (67) 

where v4a(x — y) ^ 5(x — y) as the parameter of the ultraviolet cutoff A tends to infinity. Eq. (|67D allows 
us to regularize the canonical commutation relations between field and momentum: 

bA(x),7r(x)] = Zv4a(x - y) 

instead of 

[(^(x),7r(x)] = Z(5(x-y). 

To perform infrared regularization, we consider the theory in the box with sizes L x L x ... x L with the 
periodic boundary conditions. 

After regularization eq.(p^ is written as 

z^M-)] = Jd^ I^IA^) + ^(Vv.a)^(x) + ^y^Kx)^ (68) 

where 

Mi(x) =m^ + A($i, ^l{^)^t), i^t, ^t) = 1. (69) 
Let us investigate the system of equations (|68D , (|69|) and consider the divergences in this system. 



3.4.2 Renormalization in the Schrodinger picture 



Consider the problem of renormalizing the obtained equations. For the spatially homogeneous classical 
field, analogous problem was considered in while some examples for the inhomogeneous field case 

are presented in |3^. Another approach was considered in refs. The idea is that one should 

impose the conditions not only on the counterterms entering to the Hamiltonian but also on the state 
vector. One cannot consider the state vectors which are regular as the cutoffs tends to infinity. One 
should specify how the state vector depends on the parameters of the cutoffs; this dependence is singular. 
Then one should check that these conditions are invariant under time evolutions. 

Let us perform the renormalization procedure in such a way for Gaussian and non-Gaussian solutions 
to eq. 



3.4.3 Gaussian and non-Gaussian states 

We are going to investigate eq . (|68D . First of all, suppose that the function U((x) is known. We will require 
this function to be non-singular: 

lim Mt(x) < oo (70) 

A^oo 
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because it enters to classical equations of motion for the model (|68|). 

Consider the wave functional being a product of a polynomial by the Gaussian exponent: 

J2 J cixi...cix„/„(xi,...,x„)9?(xi)...9?(x„)exp dxdyip{x)R{x,y)Lp{y) 



(71) 



To investigate the divergences in the solution of the Cauchy problem for eq.(^), it is convenient to notice 
that the operator 

1 6 



B+[Pt,Qt] = Jdx 



^(x)P*(x) - 



i 5^9 (x 



■Qt(x) 



(72) 



commutes with the operator i-^ — H {H is the operator entering to the right-hand side of eq.(^)) if 



4(x) = Pi(x), 

A(x) + I dydzAAix - y)(-Ay + utiy))AAiy - z)Qtiz) = 0. 
The initial condition ([7I| ) can be presented as a linear combination of the states like: 

B^[P^,Ql]...B^[P^,Q^]<l>, 



(73) 



(74) 



where $ is a Gaussian state. The solution to eq.(^) which obeys the initial condition (^) is expressed 
via the solutions (P^^(5^), (Pf ,Q™) to the system eq. (|7|), 



B^[P,\Ql]...B+[Pr,QT]<^K 

Since eqs.(^3]) have regular limits as A ^ oo, the operators (^) are regular as A ^ oo. This means that 
singularities may arise in the Gaussian state $* only. 
Let us consider the Gaussian solution to eq.(p3): 



Qexp 



dxdy(p{x.)Rt{yi,y)ip{y) 



(75) 



Substitute formula ( ffSD to the Schrodinger equation (|68|). By R we denote the operator with kernel 
P(x,y): 

(P/)(x)= /ciyi?(x,y)/(y). 



Analogously, by A we denote the operator with kernel A\{x, y). We find that the Gaussian quadratic form 
obeys the nonlinear equation: 

+ RtRt + i(-A + Ut)A = 0, (76) 

dt 

where Ut is the operator of multiplication by 'Ut(x). The prefactor q satisfies the equation 

|(/nQ) = -^TrPi. (77) 
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3.4.4 Operational calculus 

To investigate singularities in eg. ([7B|) , it is convenient Q to consider the symbols of the operators, 
i.e. to present each operator via operators x and —i-^, for example: 

1 

4 = i?t(x,-*^). (78) 

Since the operators x and —i-^ do not commute, it is important to specify ordering of the operators. 
Notation ([78|) means that momenta operators act first and multiplication operators act next, so that the 
operator (|78D transforms the vector 

/p(x) = conste''^'^ 

to the vector 

(4/p)(x) = /?t(x,p)/p(x) (79) 

Eq.(|79|) can be treated as a definition of the symbol of the operator Rf. 

Let us consider the product of the operators B = i?(x, — i^) and C = C(x, — i^). It can be also 
presented via coordinate and momentum operators: 

BC = {B*C){^,-i 



1 

2 . d 



It happens that the "product" B * C can be presented as 

1 

(5*C)(x,k) = 5(x,k-i^)C(x,k). (80) 

To prove formula (pOD, one can consider first the case 

C(x,k) = e*P"C(k). (81) 

In this case 

^^(x) = C(k)/k+p(x), 
5(7/k(x) = E(x, k + p)C(k)/k+p(x) = i?(x, k + p)e^P-C(k)/k(x), 

so that 

(B * C)(x, k) = 5(x, k + p)C(x, k), 
1 

5(1, k - i-^)C{x, k) = B(x, k + p)C(x, k) 
ax 



Formula (|80D is then proved for this partial case. To check eq. (|80| ) for general case, notice that if eq . (pO|) 
is valid for operators C = Ci and C* = C*2, it is also valid for C = C*i + C*2- But any function can be 
presented as a linear combination of functions (0). Thus, we justify eq.(^) for general case. 
Eq.([75D can be considered as an equation for the symbol of the operator Rt: 

i?t(x, k) + {Rt * i?t)(x, k) + * (k^ + Mt(x))Ak = 0. (82) 

We have taken into account that the operator A can be presented as A(—i-^), where A(k) = A^ is a 
Fourier transformation of the function A\. As A — oo, A^ — > 1. 
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3.4.5 Singularities of the Gaussian quadratic form 

To investigate the uhraviolet singularities, it is important to consider the behaviour of the symbol of the 
operator i?^ as k oo. It is convenient to expand eg. (pO]) into an asymptotic series: 



{B * C)(x, k) = 9^B{^^^) ^'C(x,k) _ ^gg^ 

This is an expansion in l/|k|, because the ratio between the next order of the series and the previous one 
is 0(l/|k|). We see that in the leading order of l/|k| 

(5*C)(x,k) ~ 5(x,k)C(x, k). 

Let us consider the behaviour of eq.(|82D as k ^ oo. Since the third term of the left-hand side of eq.(|82|) is 
of order |kp, one can think that Rt ~ |k| as |k| oo. Thus, it is reasonable to extract the most singular 
term from 



Rt (x, k) = z + i2^ + n (x, k) . (84) 

Eq.(P^) takes the form 



ri(x, k) + 2^AkVk2 + /iV,(x, k) + Ak(Mt(x) - = 0{\k\~'). 

We see that the quantity rt(x, k) should be of order |k|~^: 

r,(x.k)^!;l|^^^ + 0(|k|-). (85) 



Analogously, one can find a solution to the equation (|76|) with the help of the iteration procedure up to 
0(|k|~") for arbitrary a. 

Let us impose at the initial time moment the condition 

i?t(x, k) = Rr'i^, k) + Rn^, k), (86) 

where is the obtained solution to eq. (ff^) , while is of order 0(|k|""). Then eq.(|86D will be 

invariant under time evolution. One can therefore find all the singularities of the quadratic form entering 
to the Gaussian exponent. 

3.4.6 Renormalization of mass and coupling constant 



We are now going to check our main condition ([70|). It is then necessary to find the singularities of the 
average value of ip\{x.). This matrix element can be written via the functional integral 

< ^i(x) >= I D^P[^(-)]v^i(x) exp[-(^, ImR^)], (87) 

where P is a polynomial functional of ip{-) with smooth coefficient functions. Calculation of this integral 
is standard: it is presented as 



6 



Since the remaining integral is easily evaluated, the latter expression can be simplified: 



6 



/O 1 
dydzAf,{x - y)AA(x - ^) JJJJ^JJ^'^^ 
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52 



The only singularity in this average value comes from the operator 
element consists of the regular and singular parts: 

1 



This means that the matrix 



< Cpl{x) >= 2 < M|X > + < (fK^i) >reg ■ 



Since the symbol of the operator Rf is 



i?i(x, k) = zA^Jk^ + fi^ ( 1 + ^ff, + 0{\kr') 



2(k2 + 

the leading orders for the symbol of the operator A{ImR)~^A have the following form: 



{A{ImR)-'A){x,k 



+oi\kn)- 



Since the matrix element of arbitrary operator B is expressed via its symbol as 

< x|5|x >= ^ I < x|p > |25(x, p) = ^ ^ 5(x, p), 
P ^ P 

the average value (|87|) is presented up to a finite quantity < V5^(x) >reg in the following form: 



< ^li^) >-- 



1 



2Lrf ^ Vk^TT? 4Lrf ^ (k2 + /i2)3/2 
Substituting eq.(|^) to eq. (p9|) , one obtains the following condition: 



+ < ^^(x) > 



reg 



+ A 



1 



E 



4L^ Y + '"^)^^^ 



+ A < V5 (x) >reg= M^) 



Let us investigate eq. (pO]) . First of all, set 



2 2 I 



2L'^ 4- Vk2 + /x2' 



(89) 



(90) 



(91) 



Since the sum in the right-hand side of eq.(^) diverges as A — >■ cxo, one should choose the quantity 
(square of the bare mass) to be infinite to make /x^ finite. We will see in the following sections that the 
quantity yU is a physical mass of elementary particles, so that it must be finite. 
Eq.(pOD can be presented then as 



=< ip (Xj >reg 



A 



sr^ Ak 
^4(k2+/i2)3/2' 



The sum over k entering to this equation is infinite for the (3+l)-dimensional case. One should then 
choose the bare coupling constant A in such a way that 



A. 



1 



1 1 ,p 
A + L'^ ^ 4(k2 + /i2)3/2 - 



< oo. 



(92) 



The quantity Xr is a renormalized coupling constant. 
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3.5 Complex- WKB approach for the third-quantized systems 



In subsection 3.3 we have found the classical equations with the help of conjecture ( |5^ on the average 
values of the Heisenberg operators. Let us now check this conjecture and develop a systematic semiclassical 
theory. We will substitute an ansatz for the third-quantized state to the equation of motion. We show 
that the constructed state approximately satiisfies the equation. 

3.5.1 The complex- WKB ansatz 

Let us consider the following vector of the third-quantized Fock space: 

^ = exp Qs*^ U<s>tY\ (93) 

where S** is a real number, U^t is the unitary operator (0), is an e-independent state. It happens that 
the third-quantized state (p^) being a vector (^) approximately obeys the equation (^81): 



^^ = -HiV^A+, V^A-)^. (94) 
at e 



Namely, one can use the commutation relations (|5B|), as well as the formula 



^ J D^($*b(-)]A+b(-)] - $**b(-)]^-b(-)]) + yJ D^m^mnvi-)] - ^'*M-m^{-)]). 

Analogously to appendix B, one can make equal first the terms of order 0{l/e) in eq.(^): 

^ = --[($*, $*) - ($*, $*)] - H{^\ $). (95) 

Then one can consider the terms of order 0{e~^^'^) and obtain the classical equation of motion (0). Finally, 
the remaining non- vanishing as e — ^ terms of order 0(1) give us the following equation on Y^: 

= H,Y\ (96) 

where 



(97) 

An interesting feature of eq. (|97|) is that it contains terms with two creation operators. This means that 
vacuum vector is not a solution of eq.(0). Therefore, the product ( |64D is not an asymptotic solution to 
the A^-field equation: one cannot suppose the fields to be independent, it is necessary to take into account 
the correlatins between fields even in the leading order of 1/A^-approximation. 
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3.5.2 Fixing number of fields 



The constructed asymptotics ( P^D is a vector of the third-quantized Fock space. All its components are 
non-zero. However, our purpose is to construct the asymptotic solutions to the A^-field equation. This 
means that we should consider the A^-th component of the vector (^). Consider the operator P/v of 
projecting on the A^-field subspace: 

P^(^o,^i,...,^7v,...) = (0,0,...,0,^jv,0,...). 
According to subsection 3.2, the vector 

PnU^Y (98) 

is not exponentially small if 

J Dipm^{-)]\'^ = eN. (99) 

The fact that we are interested in the A^-th component of the vector (|93|) only means that expression ( p3D 
contains more information than necessary. This implies that two different vectors Yi and Y2 may lead to 
equal A^-field states PnU^Yi and P/vf/$F2- 
Namely, let us use the identity: 

Pn (^e J DipA+[^{-)]A-[ip{-)] - eN^ U^X = 0. (100) 

Applying the commutation relations (^), one obtains that 

P^U^ [^j D^\^ip{-W -eN + Vi(a+ + a') + e j D^A+[ip{-)]A-[^{-)]^ X = 0, 

where 

'Dy.A+[y,(-)]$b(-)], 



a- = j DipA~[ip{-W[v{-)]. 
Taking into account eq. ( P-OOj ), one obtains that for 

Y = {a^ + a-)X 



the X-field state (|9q) is small as e ^ 0. This means that there is an invariance of the vector ([93| ) under 
transformations 

Y ^Y + {a+ + a-)X. 

One can perform the "gauge-fixing" procedure by imposing on Y the additional condition, for example, 

a-Y = 0. (101) 

Another approach is to consider the gauge-invariant generalized state vector 

Z = 6{a+ + a-)Y (102) 

instead of the Fock space vector Y. Eq. (|102|) is associated with the vector Y by the relation 

Z = e-¥^''^Y. (103) 

Since the operator + oT commutes with the operator ij^ — H2, we find that the generalized vector 
obeys the evolution equation: 

= H,ZK (104) 

at 
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On the other hand, the vector does not obey eq. (|104|) since the operator a does not commute with 
— H2. However, the most suitable form of the A^-particle wave functional uses the vector Y . Since the 
vector Ui^Y can be expressed via components of Y as follows, 

00 -| „ _ 00 —N/2 



n=0 V n\ " j^^Q 

the expression for the A^-th component of the asymptotics is: 



N 1 

vi/*^bi(-),...,V^„(-)] = e'"^"'E^T7^ E U^n{■),■■■,V>^M n ^'Iv^i-)]- (105) 



Condition ( |1U1| ) means that the functional Y^ is orthogonal to $*: 



Thus, eq.( |105| ) has the following meaning. The n = 0-term corresponds to probability amplitude that 



all the fields are in the state $, the n-th term specifies that amplitude that N — n fields has the wave 
functional $ etc. 

3.5.3 Renormalization of the cosmological constant 

We have seen in the previous subsection that the phase factor in the solution to the classical equation 
diverges. However, the quantity also diverges. Let us show that these divergences cancel in a vector 

First of all, it is necessary to modify the classical equation of motion, eq.(^), by adding a cosmological 
constant S to the right-hand side: 

z^M-)] = jd^ Q7r^(x) + ^(V^a)^(x) + ^^i(x) + 8^ (106) 

This corresponds to usual renormalization of the vacuum energy in the quantum field theory. Modification 
( |106[ ) means that the A^-field Hamiltonian is modified by N£. Eq.(|77D for the singular pre-exponential 
factor becomes the following: 

j^{lnct) = -\TrRt - I J d-K£. (107) 
It follows from eqs.(p5|) and ( p.07|) that the quantity 

at = cte'^' 

obeys the equation: 

^{Inat) = —TrRt - i / rfx^ + ^ / rfx < (/?a(x) . (108) 



dV 2 'J 4 

Let us find the singularities of TrRf It is necessary to expand the symbol of Rt up to ©(l/lciJkl^), where 

i?i(x, k) = tA^u^ + n + r2 + rg + Oil/lu^l"), (109) 

where = O {1 / \uj]i\"^) . 

Substitute expansion ( |109| ) to eq.(^). Making use of eq.(|^), we obtain the following relations. The 
terms of orders 0(1), 0(l/|a;k|) and (^(l/lco'kP) are: 

2iA]^LJi,ri + * (m - yU^)Ak = 0, 
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d dvi 
OK ax 

Since trace of the operator is expressed via integral of its symbol over x, while integral of the derivative 
vanish, one finds: 

d Ti 
Trro = — -Tr- 



dt 2iuJk 
Trr^i = -—7:Tr — \ — Tr — —. 

We see that the singular part of TrRt is presented as: 

for some singular function depending on 'Uj(x) and its derivatives. One can then extract singularities 
from the factor at- 

at = e^bt, 

where bt is regular. The fact that the prefactor in the Schrodinger representation is singular is usual for 



quantum field theory However, this singularity vanishes if one changes the representation PB 
Thus, we see that the cosmological constant should be chosen in order to make the quantity 

1 ^ , /. i Ut(x) — i (ut(x) — A 9/ n 9 ^ 



finite. Making use of the relation (|6 



< VA^) >= + — — 

and eqs.(pID, (|92D , we obtain that the quantity —S considts of constant and regular parts: 



If one sets the bare cosmological constant to be 



the divergences in the vacuum energy will be canceled. 

3.6 Non-symmetric solutions to the large-N theory 

The formalism of second quantization works in quantum mechanics only if we require the wave function 
to be symmetric or antisymmetric. Analogously, the bosonic third quantization procedure is applicable to 
the A^-field theory only if we require the wave functional "^Nlfii'), '^n{-)] to be symmteric with respect 
to transpositions of the fields. This means that third quantization procedore allows us to construct only 
symmetric approximate solutions to the A^-field Schrodinger equation. One can notice that the functional 
( |105| ) is symmetric. 
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However, one can be interested in more general solutions to the large-A^ Schrodinger equation. To 
construct such solutions, one can consider the wave functional 

^[Xi{-),-,Xk{-),^i{-),-,^N-k{-)] 

to be symmetric with respect to transpositions of the last N — k arguments fi,..., ^N-k only. 

Let us perform now the procedure of third quantizationonly for arguments if. This means that we 
consider the extended Fock space of sets 



/ ^,[xi{-),....Xk{-)] \ 

^N-k[Xl{-),-,Xk{-),Vl{-),-,^N~k{-)]) 



V 



(110) 



The (A^ — A;)-field Hamiltonian can be presented as a sum of Hamiltonians of the fields (p, x ^-nd interacion 
term: 



where 



N-k 



^N-k 



a=l '' 



1 52 1 m? \ 

■2 J^^RJ^ + -V^.(x)V^.(x) + -^.(x)c,.(x) j + 

N-k 



Xe r 



afi=l 



is the (yj-Hamiltonian, 



^N-k 



Xe f 

+ ^ J J2 Xa(x)Xa(x)X;3(x)x;3(x) 

a, (3=1 

is the Hamiltonian of the x-field, while the interaction term is 

Xe " ^-'^ 



H"^^ =4-1 C?xX! Xa(x)XQ(x) <^a(x)(^a(x 
^ •' a=l a=l 



For the simplicity, the cutoffs are omitted. The third-quantized version of the Hamiltonian H^_j. was 
written in subsection 3.1, eg . (^Tf) . The Hamiltonian is not to be third-quantized. The interaction 
Hamiltonian is presented as 

\ k 

H^^ = ^ f rfx^xa(x)x„(x) / D^A+[^i-)]M^)M^)A-[^{-)]. 

^ •' a=l 

One can construct the asymptotics as — 00, k = const for the obtained third-quantized equation by 
the complex- WKB technique considered in the previous subsection. The ansatz (^) is substituted to the 
Schrodinger equation, where is a vector (|11CI|) . Equation for 5* and classical equation for $* coincide 
with obtained in the previous subsection. Equation for has the form: 



:ni) 
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where the operator if 2 is given by eq.(^), while the Hamihonian is 



if. 



a=l 



(112) 



We see that X" cind Fock degrees of freedom are sphtted in eg. (|1 1 1|) . Each of fields x is a free field 
interacting with the external potential Mt(x) (^) specified by the classical solution 

To fix number of fields, one can project the constructed asymptotics on the {N — /c)-field subspace. All 
further derivations are analogous to the previous subsubsection. 



3.7 Generalization to other models 

We have considered an example of the large- theory, the $'^-theory. Let us briefly discuss other large- iV 
models and applications of the third-quantized approach. 



3.7.1 The LpLfx-y^odel 

Consider the quantum field theory of scalar fields (p^,...,(p^ and one field x with the Lagrangian: 



Let us show that semiclassical third-quantized approach is applicable to this theory and find classical 
variables and classical equations for this model. 

The Hamiltonian of this model can be expressed via creation and annihilation operators [(/)(■)] in 
the third-quantized form: 



H = J DipA-^[ip{-)] J rfx 



1 5^ 1 

'2 + 2 (^^)'(^) + W + .v/^^^(x)x(x) 



1 
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When one performs the rescaling of "quantum" variables, 

A^V^=<I>^,V^X = Y, 

one finds that the Hamiltonian is proportional to 1/e, the commutator between operators $^ [</)(■)] is small, 
and the coefficient of each differentiation operator is e. This means that semiclassical methods can be 
applied to this model, while classical variables are the following: complex functional $[(/9(-)], classical field 
y(x) and canonically conjugated momentum P(x). 

The classical equations can be obtained by substituting the complex- WKB ansatz 



Uqyt exp 



I rfx[P(x)x(x) - r(x 



1 6 



to the Schrodinger equation. One will obtain the following classical dynamics: 



(113) 



c?x[ — 



Y = V, 



-V = -AY + M^r + (/($, ip{x)^). 
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3.7.2 Spontaneous symmetry breaking case 

Let us consider tlie spontaneous symmetry broken large- field theory, 

L = \d,v^d,^^ - ^(V'V'^ - Nv'f- (114) 
Since the Hamiltonian can be presented in a third-quantized form, 

= /B,A-M.)l/dx|-l^-^ 4(V,«x)lA-[,(.)] + 

+'^jd^ij DfA^[v{-)]v\y.)A-\f(.)] - J 

our semiclassical approach is applicable to this case. However, the ansatz (|93|) seems to correspond to the 
spontaneous unbroken phase of the theory because fields (p^, are not distinguished in the asymptotic 

formula. 

To specify solutions corresponding to the broken phase of the theory, it is reasonable to denote ip^ = x 
and consider the model 

One can then analyse this model analogously to the previous subsubsection; the classical variables will be 
the same: the complex functional the field Y and momentum V. 

3.7.3 0{N) - nonsymmetric theory 

The method of semiclassical third quantization can be applied to the theories of a much more wide types 
than traditional large- iV theories. For example, consider the 0(A^)-nonsymmetric $^-model 

^ = T^y.^P^d,^^ - - f - ^(E^V^)'- (116) 

Since the third-quantized form of the hamiltonian depends on the small parameter according to eq. (p^) , 
one can apply the complex- WKB approach and obtain the following classical equation: 

z^M-)] = J dA^A^) + ^(V^)^(x) + ^^\^) + g^\^))^M-)i 

where 

uti^) =m^ + <^'(x)$0, {^u ^t) = 1. 

Thus, we see that the classical theory for the large-A^ system (|116|) in the external self-consistent field. 
Although this model is not exactly solvable, one can investigate qualitative features of it and obtain an 
information about the model ( |116|) . 

4 Energy spectrum of the large-N theory 

In this section we find stationary asymptotic solutions of the large-N Schrodinger equation. 
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4.1 Approximate stationary solutions 



In the previous section we have constructed the asymptotic solution ( |1U5| ) to the A^-field equation. Let us 
find in what cases eg. (|105|) gives us a stationary state: 



(117) 



First of all, notice that the solution $* to the classical equation should depend on t as e 

= e-^^*$°. (118) 
Then, the n-th component of the third-quantized state Y is also a periodic function: 

If the conditions ( |118D and ( |119D are satisfied, eq.( |105D gives us a stationary A^- field state, because 



— = i($*, $*) - H(¥\ ¥)=n- i7($°*, <l>°), 
at 

so that the quantum energy is expressed as: 

E = NH{^*,^)+u + 0{l/N). 



(120) 



Let us investigate the obtained conditions, (|118|) and (|119| ). Condition ( |118|) means that one should 
consider the stationary analog of eq.( 



i7r^(x) + i(V^A)^(x) + ^^i(x) 



(121) 



u(x) = + A($,(^i(x)$). 
To analyze eq. ( |119| ), notice that it can be presented as: 

^—i{nQ+uj)t ^0 



(122) 



where the generalized third-quantized vector Z obeying eq. (|104D is related with Y by eq.( p.03| ). Eq.( 122 ) 
means that 

{H2Z)n = {nil + iu)Zn- 
Making use of the form of the operator of number of fields, one obtains: 



H2-nJ DipA+[ip{-)]A-[ip{-)]^ Z = ujZ. 



Eq.( |102| ) implies the following relation: 

' (A+[v.(-)]<f [V^(-)] + A-[^(-)]$*[^(-)]) ^ = 0. 



(123) 



(124) 



Thus, the problem of finding the large- spectrum is reduced to the problem of finding a solution to 
classical equation ( |121| ) and oscillator system ( |123[ ), ( |124D consisting of the oscillator Schrodinger equation 
and aonstraint relation (|124|) . 
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4.2 Ground state energy in the large-N theory 

Let us consider the following solution to eq.( |121| ). Suppose m(x) to be a constant: 

■u(x) = fj?. 

Choose the classical solution $ to be a vacuum state of the one-field system with mass ^. Thus, our main 
equation ( |121j ) takes the form: 

= m'^ + X< 0|^i(x)|0 > . 
The firld v^(x) can be expressed via creation and annihilation operators, 

^(x) = ^ E ^Ke-^'^^ + a^e^^% (125) 

where the factor ^/A^ arises because of regularization considered in eq.( |121| ), while = i/k^ + n'^. 
Therefore, one can evaluate the average value of 9?^(x): 

<0|^i(x)|0>=iji:^. (126) 

Thus, we obtain the relation (|9TD on the parameter fi playing a role of mass of elementary excitations. 
Namely, non-symmetric excitations obey Schrodinger equation with the Hamiltonian ( p.l2| ). However, this 
Hamiltonian corresponds to the free field of the mass /i. 

Let us find the leading order of the energy, NH{^*, $) expressed via average values: 

/I 1 TD^ A 

rfx[-7r2(x) + -(VyPA)^(x) + — ^a(x)] > +- < 0|^i(x)|0 (127) 

making use of eqs.( |12lD and ( |126| ), one finds: 



2N 



^-A'lfi-JZ.^E^j l+0(l) (128) 



Since the vacuum energy is ^ I^k^k^^k, eqs.( |128| ) for the ground-state energy density has the form: 

E 



One can investigate then the dependence of the vacuum energy density on the boundary conditions. 



(129) 



4.3 Structure of the ground state wave functional 



Let us find the asymptotic ground state wave functional in the large-A^ field theory. It is sufficient to 
construct a solution to oscillator-type equation ( |1 23|) which satisfies the constraint condition ( |124] ). We 
will look for this solution in a Gaussian form: 



exp 



^ /Dy,D0A+[y,(-)]i?b(-),0(-)]^+[0(-)] 



10 >, 



where |0 > is a third-quantized vacuum containing no fields. The constraint 
condition on R: 

' R[ip{-),<P{-)W[<P{-)]D<P=-^ip{-)], 



(130) 

leads to the following 
(131) 
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while Schrodinger equation ( |123| ) gives us a nonlinear relation 



+ 2 / 'i'^n (v?(x)4b.(-)l + / I'»i*?(x)*-|¥ii(-)]-Rb.(-),"?.(-)]) = 0. (132) 
It follows from eq. ( |103| ) that the third-quantized vector Y entering to eq. ( |105| ) has also the Gaussian form: 



Y = exp 



10 >, 



(133) 



where the two-field state M being equal to 

M[v.(-), 0(-)] = ^b(-), 0(-)] + <^'b(-p[0(-)] 

is orthogonal to because of eq .( 11311) . The ground state A^- field wave functional is presented ac- 

cording to eq.([L05|) as 



[N/2] 



^ivbi(-),-,¥'^(-)] = E 



1 



^ fnmin ^ Mb,,(-),¥..,(-)]-M[^.,,_,(-),¥^.J-)] n ^[^'^i-)] (134) 

1=0 l^^^ J l<ii^...^i2;<Af iy^ii-i2i 

We see that one should take into account the two-field correlations even in the leading order of the 
semiclassical approximation. 

Let us find the two-field states M and R which should obey the following condition: the norm of the 
operator with kernel M[(p{-), 0(-)] is lesser than 1. Otherwise, the vector (|133|) will not belong to the Fock 
space [p7| . 

It is convenient to consider the following one-field functionals: 



(135) 



where <l'[(y9(-)] is a vacuum state functional corresponding to the field with mass fi. The operators are 
second-quantized creation operators being linear combinations of field and momentum operators acting in 
the space of one-field states. 

We are looking for the solution to thenonlinear equation ( |132|) in the following form: 

R[M-),M-)] = -nM-mM-)]+ E Ck,k.p.p.$k.k.bi(-)]'^'pip.b2(-)] (136) 

kik2PiP2 

which automatically obeys the constraint (|131|) . 

Substitution ( p,36| ) gives us the following equation on the coefficients Ckik2pip2- 

E (^ki^ki + ^k2^k2 + ^pi^pi + ^p2^P2)Ckik2PlP2'^kik2[<^l(-)]'^PlP2b2(-)] + 

kik2PiP2 

+^ / d^li^l^i^)- < 0|^i(x)|0 >)$bi(-)] + E Ckik2k;^$kik2bi(-)] < 0|^i(x)|<l>k;k^ >] 
x[(^2,a(^)- < 0|v^i(^)|0 >)'^'[V'2(-)] + E Cp,p2p;p^$piP2b2(-)] < 0|^i(x)|$p,p, >] = (137) 
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Notice that the state vector y9^(x)|0 > contains only the vacuum and two-particle components, so that 

(^2,aW- < o|¥'i(x)|o >)$ = Yl *k,k, < $k,k.l^i(x)|o > . 

kik2 



This means that the expression (|136|) satisfies eq.( p.32| ) if 



(^ki^^kl + ^k2^^k2 + ^pil^pi + ^P2^P2)Ckik2Pip2 + 

+^/c?x <l>k,k2l<^'(x)|0 > + Yl Ckik2k;kJ, < 0|(^2(x)|$j^,^ >j 



X 



< <l>p,pjy.2(x)|0 > + Y: Cp.,p^p,p2 < 0|y.2(x)|$p,p, > = (138) 



The quantity Ckik2PiP2 is a probability amplitude that there are two particles of one type with momenta 
ki, k2 and of another type with momenta pi, p2. However, the full momentum should be equal to zero 
in the ground state. This means that Ckik2PiP2 ~ 5ki+k2+pi+P2- Let us denote by P and — P the full 
momentum of the fields, so that 

Ckik2PlP2 = '^'^^ki pi<^ki+k2+P'^pi+P2-P- 
P 

The matrix elements entering to eq . (11381) are calculable: 



^3) im^r^^in^- V '^^ '^% -»ki+k2x 



so that eq. (|138|) is simplified as: 

(^ki^^ki + ^p-ki^^p-ki + ^pi^^pi + ^P-pi^P-pi)"kipi + 



A ^ aMp'^P-p' \Mk'^P-k' „ „ 

E \ \ i^^.K + <k)(Vp; + <p;) = (139) 

k'lPi a/2^p;^p-p; J2^k;wp-k; 



Ck',)(Vp'i + "pi 

Thus, the equation on the functional R[ip{-), (f){-)] has been reduced to the equation on the function 
ttkp of two discred variables k, p at fixed P. To investigate eq.( |139|) , let us denote by the operator 
with kernel akp, by we denote the operator of multiplication by 

= (v4kCJk + Ap_ku;p-k), 
while will be the operator with the following matrix elemant: 



A VAkAp_k V-4p^p-p 

2L^ ^/2uJyU^I^ ^y'2uJpUJp_p 

Eq. ( |139[ ) takes the form: 

T^a^ + a^T^ + (1 + a^)SP(l + a^) = 0. 
The solution to this operator equation can be written as: 

= {1 + M^)-\l- M^), 
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where 

MP = (TP)-V2[(2^P)2 _^ 2(tP)1/2^P(TP)1/2]1/2(tP)-i/2_ (140) 

We see that the solution to eq.( p.32| ) is constructed if the operator entering to eg. ( P-401) is positively 
definite, 

"(rP)2 + 2(TP)i/25P(TP)i/2]V2 ^ ^ ^^^^^ 

In this case, one can construct an approximate stationary solution to the large- Schrodinger equation. 

Usually, one is able to find semiclassical eigenf unctions around the static classical solution if it is 
stable due to small perturbations. We have seen that investigation of the nonlinear equation gives us an 
alternative way to investigate this problem. We have also obtained the criteria of stability of classical 
solution (|118|), eq.(|141|). 



4.4 Perturbatively excited states 

Let us consider now the elementary excitations around the vacuum solution. It is necessary to construct 



excited eigenvectors of eq.(123). To construct such a third-quantized vectors, let us find such an operator 
A+ that 

[A+,H,-nJ DcpA+[^i.)]A-[cpi.)]] = -/3A-^ (142) 

and 

[A+, / D^iA-'[^{.m^{.)] + A-M-WM-m = 0. (143) 



This operator will satisfy then the following properties. Let Z he a solution to eq.( |123| ) which satisfies the 
condition (|124|) . Eq.( |143|) implies that the vector A~^Z will also obey the constraint (|124|) . It follows from 
eq.(p.42|) that the vector A+Z will be also an eigenstate of eq. (|123|) with eigenvalue 



Thus, the operator A"*^ shifts the energy by (3. 

Let us find such operators A"*". We are looking for them as follows: 



A+(F,G) = D^iA^M.)]G*M.)]-A-M.)]F*[^i.)]) 



(144) 



where G and F are functionals to be found. Since the Hamiltonian entering to the left-hand side of eq.( |123|) 
is quadratic, ehile operator A+ linearily depends on creation and annihilation operators, the commutator 



will be also a linear combination of these operators. Thus, eq.(142) will imply the system of two equations: 



a;'(x) 



MX, 



V.i(x)] - fi)F[^(-)] + 



A 



V- 1 rfxv,i(x) / D<p</^i{^){F[<p{.)mm+G[mm-)])m-)] = -m^i-)], 



(145) 



(/rfx[i7r2(x) + l(V^, 



X -f 



MX 



^i(x)] - ^])G[vp(■)]+ 



A 



V- I rfx^i(x) / D#i(x)(F[0(-)]$*[0(-)] +G[0(-)]$[0(-)])<^>*b(-)] =/?G[vp(-)], 



Thus, we see that spectrum of the variation system (|145|) coincides with the spectrum of differences between 
energy levels, because the operator A+ increases the energy by /3. 

An important feature is that one muct check the property A+Z ^ 0. Otherwise, the operator A+ will 
give us a trivial solution to eq.(|T2 
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When one applies the operator (^4J) to the Gaussian state ( |13CI|) , one will obtain the vector 

I DipA+[cpi-)]Y[cpi.)]Z, (146) 

where 

This functional is certainly non-zero if 

(G,G)>(F,F) (147) 

because the norm of the operator R with the kernel R[ip{-), 0(-)] is not greater than 1. Thus, the solutions 
of the system (|145|) which satisfy the condition (|147|) will give us the spectrum of elementary excitations. 

On the other hand, let us consider the solutions to eq.( |145| ) which do not obey eq.( |147[ ). Such solutions 
can be obtained if one substitute to the system ( p.45| ) F* instead of G and G* instead of F. The corre- 
sponding operator A+[G*,F*] denoted as A^[F,G] is conjugated to operator ( [144| )- Consider the vector 
A"[F, G]Z. It has the form (|4|) where 



Y[^{.)]=^F[v{.)]~^ J D0i?[^(-), 0(-)]G[0(-)]. 



However, one of the definitions of the operator R with the kernel R[ip, 4>] is the following: RG = F, see 
appendix B. This means that A^[F,G]Z = so that we should take into account only those solutions to 
the variation system that obey eq. (|147D . 

Let us find these solutions. First of all , notice that the functionals 

i^ki...k„ = 0, Gki...k„ = a+ ...a+^r, p = Ak.^ki + ... + ^k„c^k„. (148) 

obey eqs. ( P-45|) if n 7^ 0,2. Namely, the vector ip\{-x)\0 > contains vacuum and two-particle components 
only, so that the matrix element < F + G\ip\{'x)\0 > entering to eq.( p.45|) vanishes. 
One can interpret the quantity /9 as A ^ 00 

+ ...Wk„ (149) 

as an energy of the elementary excitation of the A^-field system. Why elementary? One can be surprised 
that the theory ( |149| ) is the energy of the system of n particles. However, one excitation (|148D with energy 
( |149|) is not equivalent to n excitations with energies Uk-^,..., u^^. Excitation ( |149| ) describes n particles 
of one type, while n excitations describe n particles of different types. The third-quantized language is of 
course very unusual. 

The only non-trivial elementary excitation corresponds to two particles of one type. Let us look for 
the corresponding solution to the variation system as: 

G = go^* + E ^?kik2$kik2, 

kika 

^ = /0$ + E /klk2$kik2, 

kik2 

where $kik2 is defined by formula ( |135| ). Eqs. ( |145| ) imply the following nontrivial relations: 



(AkjU^ki + ^k2l^k2 + Pj/kik2 + TTTd / '■) - / ' ■ ) =l/piP2'^Pl+P2,kl+k2 + 5'piP2"pi+P2+kl+k2j — U, 

P1P2 V^^kit^k2 V2C^P1^^P2 



, . \/^ki^k2 \/^pi^P2 X r ^_n 

(^ki^^ki + ^k2t^k2 - Pjfi'kik2 + TTTd / o . ' / O =(^PlP2'^Pl+P2.ki+k2 + /piP2'Jpi+P2+kl+k2j — U. 

P1P2 V ^^kit.^k2 V^^P1^P2 
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The next step to simplify the system is to consider the momentum of two particles to be equal to P: 



/kika — /-ki^ki+ka+PjS'kika — fl'ki^ki+ka- 

The obtained system 



(^ki^ki + ^p-ki^p-ki - /9)^ki + 7^Y1 \o ^ ^) ^ (■^Pi + 3pi) = 

2L'^ p1^2 V^Wki^P-ki V^^Pl^P-Pl 

allows us to express the functions /k and (7k via eigenvalue /3 and unknown constant ap: 



v/^ki^p- 


ki y'^pi^p- 


Pi 


V^^ki^^P- 


-ki V-^^^pi^^P- 


-Pi 


V^^ki^p- 


-ki y^^pi^p- 


-Pi 


V2cc;kiu;p- 


-ki V'^^Pi^P- 


-Pi 



if PI +9pi) = 0, 



A apV^kAp.k 1 

Jk ^ 



k 



2L'^ A/2u;kt^p-k ^kt^k + v4p_ka;p_k + P 
A apA/v4kv4p_k 1 



2L'^ A/2u;kt^p-k -^k^^k - ^p-k^^p-k + /3 
which is determined from the condition 

fp + 9p 

p V^^p^p-p 



«p = E^^#=^v/^P^p-P 



implying equation on j3: 



1 ^ ApAp.p Apu;p + Ap_pu;p_ 



A 21^^ ^ cUpCL^p-p /?2 - {ApUJp + Ap_pu;p_p) 



For d + 1 = 4, the left-hand and right-hand side of this equation diverge as A — oo. However, making us 
of the definition of the renormalized coupling constant, we obtain the regular relation: 

^ - ^ ^p + g^P-P 1 . 

Xr 2L'^^^a;pCUp_p/52-(cUp + a;p_p)2 2a;3^- ^ ' 

Eq. ( |15CI| ) allows us to evaluate the energy of bound two-particle states and scattering amplitudes for these 
particles. 



4.5 Non-perturbatively excited states 

In the previous subsection we have considered the vacuum state of the large- theory and states with finite 
number of particles. However, one can also investigate other excited states which cannot be described by 
the perturbation theory of the previous subsection. 

To understand the problem, consider the simple analogy. In ordinary quantum field theory at small 
values of the coupling constant one can apply the perturbation theory and obtain that excitations corre- 
spond to free particles of given mass in a leading order of perturbation theory. On the other hand, one can 
also apply the soliton quantization approach which is also applicable at small values of coupling constant 
and obtain non-perturbative series of states corresponding to the soliton. 

In the large- theory, the constructed asymptotic energy levels are found provided that the number 
of particles is much lesser than A^. If we consider the state consisting of A^ excitations we should take into 
account the interaction between the excitations, so that another approach is necessary. 
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To consider such a non-perturbatively excited state, it is necessary to consider another solution to the 
classical equation ( |121| ) rather than vacuum. For example, the solution 

a^\0 > 

(|0 > is the vacuum of the field system with the mass M, is a creation operator) to eq.( |121| ) will 
correspond to the quantum state of particles of a different type with the momentum k. One will also 
be able to find the corresponding Gaussian state, solutions to the variation system etc. 
Let us suppose 

m(x) = 

and 

^ = cai^a^jO> (151) 
where c is a normalizing factor. Analogously to section 4.2, one can evaluate the average value of (f\{x.): 



Eq. (121) on the mass M is to be renormalized. Making use of the definition of the physical mass of 
perturbative excitations, one obtains: 



For d = 1,2 the right-hand and left-hand sides of this equation are finite. For d=3, it is necessary to 
perform the renormalization of the coupling constant: 

M2-/i2 1 ^ 1 _ 1 M^-/.^ 1 " 

^ . /1,2 , T\/r2 . n^2 I ,,2 ~'~ 0(^^2 I ,,2^3/2'' ^ T d 2^ 



Xr 2L'^^ "Vk2 + M2 v/F + T? 2(k2+/i2)3/2^ L<^ ^ A 2 ^ 



The quantity M plays a role of a mass of small excitations around the background state corresponding to 
the classical functional (|151|) . It differs from /i. 

To find the leading order of the energy of the quantum state, let us use eq. ( |1271 ). One obtains: 



^ = n-^j d^< (x) 

so that 



E 1^. — - " . rz — - XL^fM^-m'^V 



- = -X:^kv/k^TM^ + EAk.v/kf + M2- ^ y ^ J ^ 
where we have used eq.( |121| ). The difference between this energy and ground state energy ( |129| ) is 

k i=l 

since 

The energy of non-perturbative excitation is 
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(153) 



It is interesting that for d > 5 the sum entering to eg. (|153|) is divergent. This corresponds to the fact that 
the 0^-theory in such dimensions is not renirmahzable. Usually this feature of the theory is important 
in calculations of the higher orders of perturbation theory. We see that for the method of the large- iV 
expansion the feature of the renormalizability is very important even in the leading order. 

One can also consider small excitations around the considered non-perturbative state. One should 
investigate then the variation system which has the form (|145|) where $ is not the vacuum state functional 
but the functional ( p.51| ). There will be trivial solutions to eg. ( |145| ) with F = 0, this is the case 

<G*|v5i(x)|0>=0. 

In this case the eigenvalue /3 will coincide with one of the eigenvalues of the operator Hm — ^ {Hm is the 
Hamiltonian of the field of mass M). 

To find non-trivial values of [3, denote by F and X the one-field states corresponding to the functionals 
F [(/)(■)] and G* [(/)(■)]. Egs. (p5|) will take then the form: 

(Hm -n)F + ^l rfx^i(x)(< X|^i(x)|<|. > + < $|^i(x)|F >)$ = 

(154) 

{Hm ~n)X + ^l rfx^i(x)(< $|^i(x)|X > + < F\^l{^)\^ >)$ = PX 

Let K = ki + ... + k„. To simplify the system (|154D , consider the stste X to have momentum P and F to 
have momentum 2K — P: for momentum operator P = Z^kkctkCtk '^^^ 

Px = PX, Pf = (2K - P)F. 

This condition implies that 

< X|<^i(x)|$ >=< X|e-''%i(0)e*^"|$ >= e'^^-^)" < X|<^i(0)|$ >, 

while 

< $|</^i(x)|F >= e*(K-P)x ^ $|y,2^o)|F > . 

Analogously, one can express the operator / (ixe*'*^''<y9^(x) via the projection operators IIq on the subspace 
corresponding to a given momentum Q, 

I rfxe^(^-P)>2 ^x)$ = L'^n2K-p¥''(0)$. 

Thus, one simplifies the variation system as 

{Hm -n + (3)F + __n2K-p<^i(0)<l> = 0, 

{Hm -n-/3)x + ^np^i(0)<l> = 0, 

a=<^\^l{0)\X> + <F\cpl{0)\^>, 
so that the nontrivial eigenvalues (3 are expressed from the eguation 

^ = < $|<^i(o)[(/3 - Hm + ny'iip -{(3 + Hm- fi)-^n2K-p]^i(o)|$ > 

analogously to the previous subsection. 
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5 Conclusions 



The method of second quantization is very useful in quantum many-particle mechanics even in the case of 
a fixed number of particles. For example, this approach allows us to introduce a notion of quasiparticles 
which can be created or annihihated even if all the particles of the system are stable. The quasiparticle 
conception is very important for the condensed matter theory. 

Analogously, we have seen that the notion of third quantization can be applied to the theory of 
fields and allows us to construct the approximate solutions to the A^-field functional Schrodinger equation. 



These asymptotics are expressed via the solution to eqs. (|52D and (WA)- 

Although eq. ( ^2] ) is a classical equation, it resembles a functional Schrodinger equation of quantum 
theory of one field rather than a classical field equation. Eq. (|97D is an equation for the vector of the third- 
quantized Fock space vector. The corresponding Hamiltonian is expressed via the operators which 
can be called as operators of creation and annihilation of the "quasifield" since the analogous operators in 
quantum statistics create and annihilate quasiparticles. The quantum theory of fixed number of fields is 
reduced to the theory of variable number of fields (|97|). Since the Hamiltonian is quadratic, such a model 
is exactly-solvable by the Bogoliubov-like transformation. However, the coefficients of the transformation 
are functionals expressing via the variation system ( p.45| ). 

In quantum field theory it is also assumed that the set of elementary particles is fixed: there are 
electrons, muons etc. There were no attempts to construct a theory admitting existence of the electron 
field with some probability. However, we see that the theory (^) lells us that there are probability 
amphtudes that there are no fields, that there is only one type of particles, etc. Since the theory ( p7D has 
been shown to be equivalent to the ordinary A^-field theory, the model ( pT]) obeys all the properties of the 
field theory, for example, it is relativistic invariant etc. 

The conception that particles can be created or annihilated was very important in constructing the 
quantum field theory. One can hope that the idea that number of fields can be variable will be also useful 
in developing further theories of everything. 

This work was supported by the Russian Foundation for Basic Research, project 96-01-01544. 

Appendix A. The complex- WKB semiclassical approach 

In this appendix we briefiy review the complex- WKB method which was developed in PO . 



A.l. The WKB and complex- WKB ansatz 

Semiclassical approximation is the powerful tool to construct asymptotic solutions to the quantum me- 
chanical d-dimensional equations like 



for the wave function ijj^{x) as h —>■ 0. We have seen that equations like (|155|) arise in the large-N field 



theory, while the small parameter h may be not related with the Planck constant h: in the large- A^ theory 
n= 1/N. 

The most famous semiclassical approach is the WKB-approach which allows us to construct semiclas- 
sical solutions to eq.( |155|) of the rapidly oscillating form 



where is a real function. One can easily see that the ansatz ( [15(i| ) really satisfies eq.( p.55| ), obtain the 
Hamilton- Jacobi equation for S, the expansion for ip etc. 
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However, there exists other semiclassical wave functions that also approximately obey eq.( |155| ) as 
h 0. One more example of the semiclassical ansatz to eq.( |155| ) is the complex- WKB ansatz ||20|| : 



= ^ei''ei^'^-<^V'i^) (157) 

where is a rapidly damping function of ^, while S"*, P*, do not depend on x, P* and are 

(i-dimensional vectors. The normalizing factor h'"^^^ is a corollary of the condition = 0(1). 

At fixed moment of time the wave function (|157|) is a wave packet with the width of order 0{^/h). 
The average value of the coordinate is Q*, average momentum is P*, uncertainties of the coordinate 
and momentum are of order y/K, so that the uncertainty relation 6P6Q ~ is satisfied. Note that the 
WKB-wave function has the uncertainty of the coordinate and momentum of order 0(1). Thus, the wave 
function ( p.57| ) really determines motion of the classical particle along the classical trajectory, while the 
WKB-function (|156|) does not determine a classical trajectory. 

Let us show that the wave function ( |157| ) really approximately satisfies the evolution equation ( |155| ) 
First of all, notice that extraction of the multiplier 

from the wave function 

is equivalent to shifting the differential operators: 

ox ox 

(158) 

^;,|.ei^*e^^*(^-«') = e^^*e^^'(^-«*)(^/i|- - - P\x - Q') + P'Q'). 
ot ot 

Identities ( |158| ) for operators can be justified by applying them to arbitrary test function. Eqs. (|158|) imply 
that 

{in— - - PHx - g*) + p*g* - Hix, p* - ih—))xHx) = o. (159) 

ot Ox 
Let us consider the limit h ^ oi the latter equation. At first sight, one should simply neglect all the 
terms containing h, so that one would wonder why the function S** does not depend on x. However, one 
should take into account that the wave function x^i^) is a wave packet with width of order \/h, so that 



This means that one should not only set h = in eq. (|159|) in order to obtain a leading order in h but also 



set X = g* in this equation, so that the leading semiclassical approximation gives us the following relation 
on S': 

= P'Q' - H{Q\P'), (160) 

so that 5** is the action on the trajectory of motion of the wave packet. To find next corrections, consider 
the substitution 

X - g* = ^Vh, 

so that the average number of the observable ^ is of order 0(1). Using eq.( |16(]| ), one finds that eq.(|l59D 
takes the form 

_ p^v^ _ ^^—Q* - H{Q' + v^^, P' - *^^))/*(0 = 0. (161) 
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Eq. ( |161| ) imply that (P*, Q*) should obey the Hamiltonian system: 



while the dynamical equation for /* is the oscillator Schrodinger equation: 

(z^-iJ2)/* = 0, (163) 

where 

1 d'^H 1 d'^H 1 g 1 1 d d'^H 1 g 
^ ~ T'dQdQ^ ^ 2^dQdP~idi ^ 2~ididPdP~i^ ^ ^ 

The arguments (Q*, -P*) of the classical Hamiltonian are omitted. 

Not that for ^-dependent Hamiltonian, H{Q,P) + hHi{Q,P), the classical equations (|162|) will not 
change, while the oscillator Hamiltonian ( |164D will contain the additional term, ifi(Q*,P*). 

We see that the complex- WKB approach is a good method to justify the Ehrenfest theorem. Moreover, 
one can find not only the classical trajectory of motion of the wave packet but also the evolution of the 
shape of this wave packet. Let us show how the complex- WKB approach allows us to obtain the asymptotic 
spectrum of the Hamiltonian. 

A. 2. Stationary complex- WKB solutions 

Let us find "semiclassical" stationary states with the help of the complex- WKB approach. Since the wave 
function ( p. 5 7] ) should be stationary, the wave packet should not move, so that classical coordinate and 



momentum should be time- independent, P* = const, Q* = const. The phase factor depends on t as 

5* = const - H{Q, P)t. 
The wave function ( p.57| ) depends on t as e"^^* for 

8 = H{Q,P) + he, 

if 

H2f = ef. (165) 

We see that the problem of finding approximate energy levels of an arbitrary Hamiltonian is reduced 
to the problem of finding exact spextrum of the oscillator Hamiltonian H2- We see that the oscillator 
approximation considered in [BII] is really semiclassical. 



The procedure of the diagonalization of the quadratic Hamiltonian ( |164| ) is standard One should 
consider d independent creation and annihilation operators 

(166) 

obeying the canonical commutation relations 

[A,,Ai]=S,,,[A,,AT] = 0. 

which are equivalent to 

pik)q{i)* _ q{k)p{i)* ^ ^^^^^ 
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p{k)g(i) _ q{k)p{i) ^ 0. 

The d-dimensional vectors p^^"* and q^^^ entering to creation and annihilation operators ( p.66|) are defined 
from the relations 

[H,,A^] = ±uj,Al (167) 

Since there are 2d linearily independent operators A^, one can present 2d coordinate and momenta oper- 
ators via creation andannihilation operators. Thus, eq.( |167| ) means that the operator H2 — J^k^kA^A.'j^ 



commutes with coordinate and momentum operators, therefore, it is a c-number,so that 

k 

The "frequaences" Uk are obtained from eqs.( [167| ) being equivalent to the variation system to the classical 
equations: 

icuq = » + 



One can show that the descibed procedure of diagonalization is well-defined if and only if the stationary 
solution (P, Q) of the classical equations of motion is stable. 

The wave function of the ground state is defined from the condition 

Ar/ = o,...,A^/ = o. 

The function / has the Gaussian form 

1 

f = exp(-^aO, 

where 

p^*^) = ag(^) (168) 

one can show that the matrixconsisting of d vectors reversible, so that eq. ( |168| ) determines 

the matrix a in a unique fashion. The ground state energy can be written as 

The excited states are found by applying creation operators to the ground state 

{Atr-iAjT'^U (170) 

and have energies 

e^+^^fcn,. (171) 

We have shown that the stationary case of the complex- WKB approach leads to the oscillator approxima- 
tion. 

Appendix B. The complex- WKB method in Fock space 

In this appendix we develop the complex- WKB method for the second-quantized equations. For more 
details, see ^ . 
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B.l The complex- WKB ansatz 

Consider the Schrodinger equation 

i— = H^' (172) 
at 

for the time-dependent vector \1/ of the Fock space. It happens that semiclassical methods are apphcable 
to eq.( |172|) if the Hamiltonian depends on the bosonic creation and annihilation operators and small 
parameter e as follows: 

H = ^H{^/ea+, v^a") (173) 

For the case of a finite number of degrees of freedom one can consider the following representation for the 
creations and annihilation operators. The vacuum state corresponds to the wave function 

^o(0 = exp(-^E^')' (174) 
while the operators can be presented as 



= ' -/-^ (175) 



Substituting eq. (|175|) to eq. ( |173| ), one finds that the Schrodinger equation (|172|) is taken to the semiclassical 
form: 

dip ^^,x — d/dx x + d/dx. , 

while the analog of the Planck constant is e, e = h. 

We see that one can apply the complex- WKB method to the Hamiltonian of the type (|1 73|) . Note that 
A^-field Hamiltonians considered in this paper depend on the small parameter analogously to ( |173| ). 

To construct an analog of the wave packet ( |157|) , let us consider the partial case of the formula ( [L57| ), 
the particle eith zero coordinate and momentum, with the wave function 

/(x/v/i). 

It is easy to see that this quantum state corresponds to an e-independent state vector of the Fock space. 
Really, such a vector can be expressed via creation operators and vacuum state ( |174| ). However, the 
wave function (|174|) depends on x/y^ only, while the creation operators are expressed via operators 
of multiplication by x/y^ and operator of differentiation with respect to Xij^fe. This means that e- 
independent state describes a wave packet with P = 0, Q = 0. 

To construct the Fock-space analog of the wave function ( [157| ) in general case, present it as 

Since operators x and djdx are linear combinations of creation and annihilation operators, the exponent 
can be presented as 

for some numbers ^i. Thus, we see that the complex- WKB ansatz in the Fock space should depend on e 
as follows: 

^* = V^t^s^Y\ (176) 

where 
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This ansatz can be considered also in the infinite-dimensional case. 

Analogously to appendix A, one can show that vector ( |176| ) approximately satisfies eq.( |172| ). Namely, 
the following relations take place: 



<fk_ 



(177) 



:i78) 



t 



In the leading order 0(1/6:) we obtain the condition on S 



S='^U^l^-^-^^^)-H{^*,^)- (179) 



The next order 0{l/^/e) leads us to classical equations for classical variables ifl: 

. difk dH 
dt difl ' 

while the terms of order 0(1) allow us to obtain the equation for the Fock space vector F*: 

dY' 
i—r- 
dt 

where 



:i8o) 



H2Y\ (181) 



1 d'H , , d'H _ 1 _ d'H 
B.2 Asymptotic spectrum of the A^-particle Hamiltonian 

Let us consider the Hamiltonians ( |173| ) that consists of terms with equal number of creation and annihi- 
lation operators, i.e. the Hamiltonians that are invariant under the global transformations, — > a^e^*". 
This means that the numbar of particles is an integarl of motion. 

To find eigenstates of the iV-particle component of the Hamiltonian, consider the iV-th component 
of eq.( |17(j| ) and find find whether it is a stationary vector. For the simplicity, set N = 1/e. Eqs.( |177] ) 
implies that the average value of any observable /(^ Z^m c^m'^m)? where / is an e-independent function, is 
the following for the state \E'*, 

m m 

This means that the average value of eN is VmVm- To make sure that the A^-th component is not 
exponentially small, one should set 

E<^- = 1- (183) 

m 

Let the classical solution to eq. (|180|) be 

= ^ke-'""', (184) 
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while the n-th component of depend on t as 

= K^e"^"^*-^"^*. (185) 

It follows from the relation 
that in this case 

Eg. ( |185|) means that the vector e^'^'e*^*^™"™"™!^* is time-independent. Therefore, the A^-th component of 
the state t/<^t,ol^* behaves as e-^(^^+'^)*. Eg .(11791) implies S = {Q-H{<^*, (p))t, so that the A^-th component 
of the vector ( |176| ) depends on t as 

and gives us an eigenvector of the A^-particle Hamiltonian with the energy 

E = NH{ip*,ip)+uj 

up to 0{l/N). 

The next observation is that if one chooses the vector Y as 



the iV-th component of the expression ( |176| ) will vanish. This is a straightforward corollary of canonical 
commutation relations ( [177| ) and the condition ( |183| ). Since the vector Y is defined up to 0{^/e), one can 
conclude that if the more weak condition 

Y' = e-^™E™«^-^--*r + J2{a+iPm + a-OX* (187) 

m 

for some time- dependent vector X* is satisfied instead of eg.( |185| ), the X-particle wave function will be 
also stationary, so that eg. ( |186| ) will be valid up to 0{^/e). 

In order to avoid appearence of the Fock space vector X*, it is convenient to consider the generalized 
state vector 

^* = ^(E(«m^- + «mO))^*- (188) 



The vector is defined from the relation ( |188|) not uniguely but up to the state vector of the form 
J2m{(^m'^rn + cim^m)-^^- Howcvcr, different vectors F* obeying eg. (|188|) determine the same X-particle 
state PjsiU^Y up to 0{y/e). This means that the vector specifies the asymptotic formula for the 
iV-particle state uniguely. 

Eg.( |187D imphes that the vector Z should obey the following stationary eguation 

H2Z = uZ (189) 

and the constraint eguation: 

{j:i(^i^n. + a;;,ip*J)z' = o. (190) 

m 

Thus, one should investigate the Schrodinger eguation with the guadratic Hamiltonian for the constrained 
system. Analogously to subsection A. 2, introduce new creation and annihilation operators being linear 
combinations of the operators a^: 



I ^ 



I 
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(191) 



Such operators obey eq.( |167| ) if the variation system is satisfied: 



(192) 



The ground state vector of eq.( |189|) is to be found from the conditions: 

A-,Z = 0. (193) 

Eqs.( |193D and ( |190| ) determine the state uniquely if the set of vectors Lp*, G^^\ G^'^\... is complete. The 
new vacuum vector will be of the Gaussian type 

Z = exp ^ a+Ri^aA |0 >, (194) 

where the operator R entering to the quadratic form is defined from the relations: 

The fact the the set G^^\ G^^),... may be not complete is important. This means that the non-stationary 
analog of eq.( [L92| ) may have one linearily growing at t solution. If we considered the more strong relation 
on F, eq.( |185| ), we would be unable to construct an asymptotics in such a case. 
One of the solutions to eq.( |188| ) on the vector Y is the following: 



y = exp (i^a+(i?;„ + ¥?z(/?^)a+ j |0 

\ Irn I 



> . (195) 



One should then ensure that the expression (|195|) really determines the Fock space vector, i.e. that |^| 
J2mn l-RmnP < oo and ||M|| < 1, whcre the matrix of the operator M is Mmn = Rmn + ^m^n- Since 
MG^^^ = —F^''\ one should require that 

||G«||>||F(^)||, (196) 

this rule allows us to select one of two frequences, u or -u. 

Excited states are expressed by eq. (|17CI|) , while their energies have the form (|171D . We see that the 
complex- WKB approach allows us to find asymptotic spectrum of the A^-particle Hamiltonian which 
corresponds to periodic solutions (|184|) of eq. ( |180| ). For more details, see |]39| 
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